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It is useful to describe a leading order parton shower as the solution of a linear equation that 
specifies how the state of the partons evolves. This description involves an essential approx- 
imation of a strong ordering of virtualities as the shower progresses from a hard interaction 
to softer interactions. If this is to be the only approximation, then the partons should carry 
color and spin and quantum interference graphs should be included. We explain how the 
evolution equation for this kind of a shower can be formulated. We discuss briefiy our efforts 
to implement this evolution equation numerically. 



1 Introduction 

The evolution of a parton shower can be understood as a numerical solution of a hnear evolution 
equation of the forns-^ 

\p{t)) =U{t,to)\p{to)) (1) 
with ^ 

U{t, t') = M{t, t') + j^^ dr U{t, t) Hiir) M{t, t') . (2) 

Here \p{t)) describes the system at "shower time" t, where increasing t denotes increasingly soft 
interactions. The state at time t is related to the state at an earlier time to by a linear evolution 
operator U(t, to). The evolution equation ([2]) is written using an operator Af{t, t') that represents 
the probability for the system to go from time t' to time t with no splitting and another operator 



'Hi{t) that represents a parton spUtting to two partons. Thus the system either goes from to to 
t with no sphtting or else it goes from to to an intermediate time r with no splitting, then splits 
at time r, then evolves from r to t with the full evolution operator, possibly involving further 
splittings. 

In the simplest sort of shower, each parton carries a momentum p, so that m partons carry 
momenta {pi,p2, ■ ■ ■ ,Pm} = {p}m- We can denote the state in which there m partons with 
momenta {p}m by |{p}m)- Then a general state |p) is a linear combination of basis states 
Kpjm), with ({p}m|/o) representing the probability for there to be m partons with momenta 
{p}m- When the splitting operator Hi{t) acts on an m-parton state |{p}m)i it produces an 
m + 1 parton state with a definite probability 

{{p}m+i\'Hi{T)\{p}m) . 

There is a certain amount of freedom in specifying 7^i(r), but one is constrained by the structure 
of the underlying quantum theory in the limit that the virtuality of the new pair of daughter 
partons approaches zero. In this simplest sort of shower, the basis states | {p}m) are eigenstates 
of the no-splitting operator 7V(t,t'), 

A/'(t, t') I {p}m) = A(t, t'; {p}ra) \ {p}m) ■ (3) 

The eigenvalue A(t,t'; {p}m) is the Sudakov factor 

A(t,to;M^) = exp(^-^^ dr j [d{pU+i]{{p}m+i\Hi{T)\{pU)^ . (4) 

The integrand 

(^^\iy^ j [d{p}m+i]{{p}m.+i\ni{T)\{p}^) 

is the total probability for the given state to split at time r. Then A(t, to; {p}m) is the probability 
for the state not to split between times to and t. 



2 The parton shower in quantum chromo dynamics 

So far, I have described a fairlygeneral structure for a parton shower. This description might 
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apply, for example, to Pythia'^J The style of the description sketched above emphasizes that 
we are using a definite linear evolution equation, so that we can bring the full power of linear 
algebra to bear on the problem as needed. At this point, we need to expand the description so 
that it can encompasses quantum interference, spin, and color in quantum chromodynamicsP 

To include quantum interference, the description should be based on the quantum amplitude. 
The quantum amplitude depends on the spins and colors of the partons, so we start with 



M{{pju)Z'Xfs\: 



Here, for hadron-hadron scattering, the partons carry labels a, b, 1, • • • ,m and each parton has 
a momentum p, a flavor /, a spin index s and a color index c. The array M can be thought of 
as a vector in spin and color space, 

\M{{pJU)) . 

The cross section for a possibly spin and color dependent observable F, including the proper 
factors for the parton distribution functions and for the number of color states iic of each parton, 
is 



4nc(a)nc(6) 2rya??bPA-PB (5) 
X {M{{p,fU)\F{{p,f}m)\M{{p,fU)) . 



We rewrite cr[F] in the form of a trace over the color® spin space, 



^[^] = / [^{^^' M Tr{p({P, f}m)Fi{p, fU)} , (6) 

where 

Thus p is the density operator in color spin space. We can expand p{{p, f}m) in basis states 
|{s,c}rrt) for the color (g) spin space, 

/}m) = |{s,c}r„)p({p,/,s',c',S,c}m) ({s',c'}rri| . (8) 



s,c s'.c' 



Thus p{{p, f, s' ,c' , s,c}m) is the density matrix. It is this matrix, with variable numbers of 
partons m, that is the basic object that evolves in a parton shower. In analogy with the 
notation used in the introduction, we consider /9 to be a vector with 

P{{P, /, s',C, S, C}m) = {{p, f, s', c', S, C}m\p) . (9) 

Notice that each parton is described by its momentum, its flavor, two spin indices, and two color 
indices. 

With this formulation, we can defindU a splitting operator Tiiit) based on the behavior of 
the amplitude when two partons become collinear or one becomes soft. This gives a shower 
evolution equation of the form Q. However, in general the no-splitting operator AA(t, t') is now 
a matrix in the color space. 



3 Issues of implementation 

A parton shower program like Pythit^ starts with a state \p{toj) with just a few partons and 
produces states \{p, f}m) with many partons at a final shower time t{. A parton shower program 
could also report a weight w for the state. The weight times the probability to produce state 
{p, f}m is 

{{p,f}mp{ti,to)\p{to)) . (10) 

How can the evolution equation discussed above be implemented as a computer program? The 
evolution equation, solved iteratively, produces results expressed as integrals, so one could simply 
perform the integrations by numerical Monte Carlo integration, producing events and accompa- 
nying weights. However, for a large number of splittings it is likely that the fluctuations in the 
weights are too large for this most straightforward method to be practical. 

To proceed, we need to make some further approximations, with the understanding that 
any approximations should allow one to approach the exact solution of the evolution equation 
by using a sequence of approximations that become more and more exact as one proceeds, 
presumably at the cost of requiring more and more computer power. 

The base approximation is to average over spins and take the leading color approximation, 
— > 0, where A'c = 3 is the number of colors. With these approximations, we find"^ that 
the evolution equation has the proper form to be implemented as a Markov process. In terms 
of numerical integration, this means that the integrals are nested and one can take the weights 
to be 1. 

Next, we need to put spin back^l We assume that the final measurement function does not 
measure parton spin. However, the spin states of intermediate partons can influence the angular 



distributions of splittings. If we use the spin-averaged shower to generate events, then the prob- 
abihty to generate a given shower history is wrong by the ratio of the spUtting probabihties with 
spin to those without spin. We can take that ratio to be a weight that accompanies the event. 
Following an insig ht of Collinspl we find that the spin weight factor can be calculated efficiently, 
using computational resources that are linear in the number of partons. (Herwig incorporates 
some spin effects using a related methodP) Possibly, for reasons of numerical convergence, one 
should include the spin exactly for the first splittings, then average over spins for further 
splittings. Then the exact result is approximated more and more closely as we take N to be 
large. 

Finally, we need to put color back. This is more complicated than putting spin back. 
We expect to use a base approximation that is much less restrictive than the leading color 
approximation but that still allows efficient computation. There is a difference between the 
exact Til and the approximate one. This difference would be included perturbatively at whatever 
order is needed to obtain an accurate result. Our work on color is in progress. 
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